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GUIDANCE EQUATIONS FOR QUASI-OPTIMUM SPACE MANEUVERS 


By D. J. Jezewski 


SUMMARY 


Approximate time- opt imum guidance equations for planetary ascent and de- 
scent trajectories at constant thrust are derived for a flat-hody two- 
dimensional case. No aerodynamic forces are considered. In addition, this 
case was extended hy introducing power series to simulate the deleted terms 
of the round-body solution. A power series in time was chosen to retain the 
form of the guidance law. The resulting transcendental equations are not ex- 
plicit in the control variables; hence an iterative scheme is used to solve 
the conditioned equations for the required control. 

Trajectories resulting from this set of equations are compared with the 
strict optimum for a number of cases to determine the proximity of the solu- 
tions . The results match well for both state and control variables in the 
cases considered. 

The guidance application was investigated by using the analytic solution 
to steer an exact integrated model. A comparison of these results with an 
optimum solution for ascent and descent trajectories shows good agreement. 

To demonstrate its utility as a trajectory scanner, the technique is ap- 
plied to the investigation of lunar landing from an elliptic transfer orbit . 


INTRODUCTION 


The requirement for optimum trajectories in preliminary design work has 
steadily increased in the past few years . The chief deterrent to the solution 
of the nonlinear differential equations has been the lack of a process which 
would converge rapidly to the desired boundary conditions. The first-order 
perturbation techniques used in convergence processes normally work well in 
the vicinity of the solution where the linearity assumption holds true, but 
they are usually unpredictable in areas not near the solution. These, unfor- 
tunately, are the areas most frequently encountered. 

An analytical scheme which provides reasonable approximations to the 
optimum solution would be of significant value, first, as a scanning technique 
for providing rapid estimates of optimum performance and, second, as a method 
of determining approximate initial values for the true optimum solution. The 
use of such a technique as a guidance scheme also has possibilities if the so- 
lution time is not unreasonable and if the errors generated by the approxima- 
tions to the true dynamic equations are not excessive. 



The present paper deals with an analytical technique derived by operating 
on a set of linearized two-body equations with the calculus of variations to 
determine a guidance law. The planar case, in which aerodynamic forces are 
neglected, is treated. The solution differs from that derived in reference 1, 
in that for the present paper a constant thrust rather than a constant accel- 
eration is assumed. The equations are first reduced to a flat-body set by 
allowing the radius to approach infinity. By integration, three simultaneous 
transcendental equations in three unknowns can be derived from this set of 
equations . Trajectories resulting from this set are compared with optimum 
trajectories to determine what errors are present and which terms, deleted by 
the flat-body approximation, are the principal contributors to the error. 
Suitable approximate expressions are introduced to the flat-body set of equa- 
tions to simulate the deleted terms in an attempt to reduce these errors . The 
equations are again integrated, and the results are compared with the optimum 
for a number of trajectories. The scheme is tested as a guidance technique by 
using the analytic solution as a feedback to the exact integrated model to de- 
termine the speed of solution and the penalty incurred as compared with the 
optimum. 


SYMBOLS 


A constant defined by equation (20) 

a n ,a 0 , ...,a coefficients of time series 
1' 2' ’ n 


a. . elements of partial matrix 


u 


C 


l' C 2'C 3 , c 4 


components of acceleration defined by equations (24) and ( 25 ), 
ft/sec^ 

initial values of Lagrange multipliers, radians/sec 


c 


F 1' F 2' F 3 


f 

G 


g 

H 


kqj ^ 2 ’ ' ' • > k 


n 


effective exhaust velocity, ft /sec 
transcendental functions 

function defined by equation (Al) 
function defined by equation (A2) 
gravitational acceleration, ft /sec* 
modified Hamiltonian 
coefficients of time series 
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I 


performance index 


I 

sp 

K 

k 

M 

m 

P 

Q 

q(t) 

r 

r b 

T 

t 

u,v 

V 

V 

c 

W 

y s 

z 

V a 2 

r 

At 

AV 


specific impulse, sec 

constant defined by equations (22) and (23) 

function defined by equation (21) 

boundary vector constraint 

mass of vehicle, slugs 

functions of time 

functions of time 

state vector 

radius to vehicle from center of body, ft 
radius of the reference body, ft 

thrust, lb 
time, sec 

components of total velocity in x and y directions, 
respectively, ft /sec 

total velocity, ft / sec 

characteristic velocity, ft /sec 

weight of vehicle, lb 

rotating coordinates of mass particle, ft 
distance from origin to surface of body, ft 

tangent of 0 

functions defined by equations (19) and (l8), respectively 

flight-path angle, deg 
increment of time, sec 
change in total velocity, ft /sec 



AV 

c 

change in characteristic velocity, ft/sec 


e t 

tolerance on time-to-go 



tolerance on initial pitch angle 



tolerance on final pitch angle 


T| 

true anomaly 


0 

pitch angle, deg, or thrust angle with respect to x-axis, 
radians 


r(t) 

control vector 



Lagrange multiplier, radians/sec 


V* 

normalized mass ratio, ny'm^ 


V 

constant Lagrange multiplier 


&2* a 3 

functions defined by equations (C24), (C25), and (C2 6), 
respectively 
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0 

initial thrust acceleration, T/m Q , ft/sec 


*1>*3 

functions defined by equations (3 6 ) and (C27), respectively 


Subscripts: 



0,1 

initial and final values of state and control variables 


A 

analytic quantity 


I 

integrated quantity 


ij 0 

partial derivative indices 


m 

order of correction 


min 

minimum 


n 

nth order term of time series 
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Superscript : 
T 


transpose 


Operators : 

( * ) time dif f erentation 

5( ) variational operator 


FQEMJIATION OF THE PROBLEM 


The mathematical model employed in this investigation is a mass particle 
with two degrees of freedom referred to a set of rotating coordinates with 
the x-axis alined along the local horizon. The axis system and the associated 
notation are illustrated in figure 1. 

The equations of motion for this model are 


where 



(l) 


(2) 

cos 9 - uv/r 

(5) 

2 / 

sm 0 - g + u /r 

W 

= 1 + p.t 

(5) 

II 

B • 

<?* 

(6) 

" s o( r t > /r ) 2 

(7) 

II 

o' 4 

1 

CQ 

(8) 
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The vehicle is assumed to have a constant thrust and mass flow rate. The 
angle 0 is the control variable. 

The final time t^ was selected to be minimized since under the pre- 
viously mentioned assumptions it will yield minimum fuel consumption. 

If in equations (3), ( 4 ), and (7), r^ is allowed to approach infinity, 
the coupling effects vanish and the gravitational acceleration approaches its 



surface value g Q . The following flat body equations result: 

X = u (9) 

y = v (10) 

u = ^ cos 9 (ll) 

V> 

v = sin 9 - g 0 (12) 

The following pitch-angle relationship (ref, l) is obtained when the 
calculus of variations is applied to this set of equations and the final range 
is allowed to be free. 


tan Q = 


C 4 ■ C 2 t 


(13) 


The derivation of the guidance law (eq. (13)) is given in appendix A. 


The solution is now entirely specified, because when the constants 
C 2 /C 5 and C^/C^ are given, the control variable is specified; and 

equations ( 9 ) to (12) may be integrated to yield a solution. To determine 
a unique solution, the constants C must be found as functions of the de- 
sired boundary conditions. This result is accomplished by integrating 
equations ( 9 ) to (12) between the bounds specified at t = 0 and t = t 
The constants C^/C^ and C^/c^ are otrtaineci ^ evaluating equa- 
tion ( 13 ) at the boundaries and yield 


1 * 


cj = taja e o 


tan 0 q - ta n 9-]_ 


( 1*0 




where 0(0) = 0 n and 


e(t 1 ) = 9 r 


Integration of equations ( 9 ) to (12) results in the following three sim- 


ultaneous transcendental equations in the three unknowns 


9 O’ 9 1’ 


and t 


1 * 


7 


I 



/ (xt-A a 2 / 

l i v i - + —) ' T v 


cp sec 0^ - A ) =0 


(15) 


, s o\ ^i K 

y l - >*l y 0 + IT " 1 " -2 


log 


M* L 


sec 0 n \ Kjj,a, 


/tan 0^ + K sec 0^ 
el tan 0~ + K sec 0 n 


+ Kttjsec 0 Q - j - - “ 0 ) 


= 0 


(16) 


Kk 1 ^ 2 (J- 

9^ 


1 ( U 1 - u 0 J - log e 


k - 


Vl 


l/2 

(Xg Kk ' sec 0^ 


1/2 ~ 

k - - a 2 Kk 7 sec 0 ( 


)H 


= 0 


(17) 


The integrations and associated transformations are indicated in appendix B 
In equations (15) to (17); 

|i.t 

crt = 


"2 ^tan e 1 - tan 0 Q ^ 


ct^ = 1 - a 2 tan 0 Q 


(18) 

(19) 


A = - Sr ( v 0 - ^ y o) + 9 sec e o 


( 20 ) 


k = a^ 2 + a 2 2 (2l) 

K = 1 (“ f ^ 6 ^ f) ( 22 ) 

K = -1 (g< 0 < (23) 

Because the control quantities t^ 0 Q , and 0 1 occur implicitly in this 

set of equations, it is impossible to solve for them directly. Since the 
equations are transcendental, a first-order perturbation technique (ref. 2) is 
utilized to force the solution to converge. This technique, which generates a 
partial mat rix, assumes that a solution exists if the matrix of partials has 
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J 



an inverse. This convergence technique and the associated partial derivative 
matrix are given in appendix C. 

Figures 2(a) and 2(h) illustrate the errors of this solution as compared 
with an optimum integrated round body for different values of the constant g^. 

It should be noted that the errors generally decrease with increasing T/Wq. 
This follows logically since burning time t^ oc (t/w q)" 1 and surface inte- 
grated over approaches a flat body as the burning time approaches zero. 

The question of what can be done to improve this solution (indicated in 
figs. 2(a) and 2(b)) without destroying the analytic properties of the problem 
is now posed. Certainly, the errors are the result of neglecting the cross- 
coupling effects and setting g equal to a constant . 

NASA-S-64-285 



(a) Characteristic velocity error. 


Figure 2.- Flat-body error for a launch to lunar orbit. - 420 sec. 
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(b) Initial pitch angle error. 

Figure 2.- Concluded. 

Figures 3(a) and 3(b) illustrate time histories of the terms 

a^ = uv/r (24) 

a v = “ g + u 2 /r (25) 

* 

for an optimum launch to lunar orbit as compared with those of the flat-body 
approximations. The displacement between the curves in figures 3(a) and 3(b) 
represents the horizontal and normal acceleration errors of the flat-body 

2 

•approximation and amounts to a maximum of 5*5 ft /sec for the normal component 
at the terminal point. The horizontal acceleration term a^ for opt imu m time 
descents or ascents from orbit can be considered to have only a second-order 
effect. For the case illustrated, a launch to lunar orbit, the max imu m value 
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Horizontal acceleration term, a u# ft/sec^ 
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(a) Horizontal acceleration tern. 


Figure 3.- Comparison of solution obtained by using analytic and exact solu- 
tions for launch to lunar orbit, t/w^ = 0.6; I = 315*0 sec. 

ii ' I C3 *0 


Normal acceleration term, a w/ ft/sec 


- 6.1 


- 5.0 


CN 


- 4.0 



(b) Normal acceleration term. 

Figure 3 »- Concluded. 

2 

of this term amounts to less than 0.2 ft/sec . Therefore, the quantity a 
accounts for the largest portion of the observed error. 


To offset the errors, functions are defined in the form 



( 27 ) 


l 



n=0 


which have the same values as the functions defined in equations (24) and ( 25 ) 

at the two bounds. The introduction of a function a = a (a), where q is 

y V 

any of the dependent variables, would necessarily complicate (if not completely 
destroy) the guidance law stated in equation (13) . A time series is accept- 
able because it will not alter the guidance law but will change the trans- 
versality condition of the calculus of variations. 


Based on this analysis, equations (l) to (4) are 


x = u 


y = v 


u = ^ cos 0 - ) h. t 

\x L_ j n +l 

n=0 


n 


v . ^ si” 8 + ^ a n+1 t 
0 


n 


where, from equation ( 13 )^ 


sin 0 = 


Kz 


"A 




K 


COS 0 7 b\ i In f 

(l + z 


5 ) 1/2 


tan 0 = z 


J 


( 28 ) 

(29) 

(30) 


(31) 


(32) 


Integrating these equations as outlined in appendix B results in the 
following three simultaneous transcendental equations. The method of solution 
of these equations is outlined in appendix C. 


13 



I , n+1 

a . t. 


»*1 


V p+ii 

r l “ ^ y 0 / , n + 1 


n=0 


X + 


jlt^n+l) 

n + 2 


a.,- 


_ — (cp sec 0 1 - A) =o 
V> 


(33) 


Kk 1 ^ |1 

9 a o 


( U 1 " u 0 + *l) " l0g e 


Vl " Ka 2 k ' 


1/2 


sec 0. 


(k - - Ka 2 k 1 / 2 sec 9 Q ^ 1 


= 0 


(3^) 


l . n+2 

— a 

n+1 l 


|J 'i y 0 y l + (n+1) (n+2) 

n=0 

/ sec 0, ' 


|iAa 2 t 1 + Wi k 


log 


tan 0^ + K sec 0^' 
e y tan 0^ + K sec 0^ 


+ Ka 0 sec 0. - 

2 \ 0 p, 1 


KP-a 

— (u - U-. + T]f. 

cp a 2 \ 1 0 1 


“I 


= 0 


(35) 


where 


% = 


l , , n+1 

h n+l t l 
n+1 


(36) 


n=0 


The required number and the evaluation of coefficients to simulate suffi- 
ciently the optimum round body trajectories must be determined. Obviously, 
the problem is reduced to a flat body if only the n = 0 term is retained 
and the coefficient hq is chosen to be 0. Consider a case where the con- 
trol program for the integrated and the analytic solutions are identical in 
time. As higher order terms in n are introduced to the flat-body solution, 
the instantaneous accelerations and hence the state variables of the analytic 
model would necessarily converge to those of the integrated solution. If the 
control program of the analytic model (eq. 13) is considered an off-nominal 
solution, the above qualitative argument implies the following generalization. 
The degree of success in utilizing this analytic technique to simulate optimum 
trajectories depends on the success in simulating the deleted accelerations of 
the flat body with power series in time. 

For a case of launch to lunar orbit (figs. 3(a) and 3(b)), it is apparent 
that the polynomials must be at a minimum cubics and certainly of higher order 
if greater accuracy is desired. This minimum solution necessitates four co- 
efficients. The results of using this truncated series to simulate the deleted 
accelerations of the round body are indicated in figures 3(a) and 3(b) . As 
was stated above, two values of a^ and (eqs. (2^) and (25)) are known. 

by their evaluation at t = 0 and t = t^ Two more solutions are required to 

evaluate the four coefficients of each series . It is apparent that once the 

trajectory progresses beyond a few integration steps in a guidance problem, 

l4 



the past history of the dependent variables will yield the required two addi- 
tional solutions. At the initiation of the trajectory, the following technique 
is used to evaluate the unknown coefficients . The derivatives of equa- 
tions (24) and (25) with respect to time are 


a 

u 


[v(u - aj + uv] 


a 

v 



> 

J 


(37) 


These equations and equations (l) to (8) are integrated either backward or 
forward in time twice by a small increment At to define the additional co- 
efficient in equations (26) and (27). These two solutions will contain an 
error, the magnitude being proportional to the uncertainty in the pitch angle 
0 at t = 0. This error in the coefficients is reduced by cycling through 
the algebraic equations and integrations until the coefficients converge to a 
set of values consistent with the computed value of the pitch angle. 


RESULTS AMD DISCUSSION 


Examples of results using the equations derived in this report are shown 
in an attempt to indicate the broad usage to which this analysis applies. An 
IBM type 1620 electronic data processing machine was used for the computations; 
the solution time for a typical boundary- value problem, as a launch to 'orbit 
or landing, was approximately 2 minutes and required five or six iterations to 
obtain five-decimal- place accuracy. Guidance problems, of course, ran much 
longer on the machine since they consisted of a series of boundary- value prob- 
lems, and convergence of solution was slow when "time-to-go" became small. 
Therefore, an open-loop system had to be used for the last portion of the tra- 
jectory. 


Guidance Trajectories 

The accuracy of guided solutions in simulating the strict optimum is 
indicated in figures 4 to 6, and terminal values are presented in tables I to 
III. In these solutions, the true dynamic equations were integrated at a 
step interval of 10 seconds, and the analytic solution was used as a feedback 
control . 

Figures 4(a) to 4(d) are time-history comparisons of the state and control 
variables for the guided analytic solution with the optimum for a launch to 
lunar orbit. Table I indicates the numerical differences at the boundary. 

This trajectory required the use of an open- loop system at a time-to-go of 
approximately 9 seconds . The largest apparent deviation in these curves is 
in the pitch angle, where there appears to be as much differential area above 
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as below the optimum. This apparent canceling is the probable cause of the 
closeness of the final values of time indicated in table I. It can be ob- 
served by examining these curves that the analytic technique has a linearizing 
effect on the solution. This appears to be a general characteristic of the 
analytic technique. 


TABLE I.- LAUNCH TO LUNAR ORBIT 

[t/N 0 =0.6; I gp = 315 sec; y Q = 1,000 ft ; v Q = 100 ft /sec; r Q = 90°] 



Optimum, integrated 

Analytic, guided 

t^, sec 

227.92 

228.87 

CD 

O 

V. 

CD 

• 

• 

• 

• 

• 

• 

27.68 

26.68 

deg 

-8.25 

-8.77 

V , ft/sec . . . . 

5,775.3 

5,807.9 

y ± , ft 

49,958 

49,962 

v^, ft /sec . . . . 

5,606.7 

5,601.9 

q, deg 

0.00084 

0.00527 

ft 

564,573- 

567,082 


Figures 5(&) to 5(&) are the time-history comparisons for a guided landing 
from lunar orbit. Numerical results are presented in table II. Optimum retro- 
grade solutions are extremely sensitive to initial conditions when solved by 
integrating the dynamic and Euler- Lagrange equations forward in time. There- 
fore, the landing from lunar orbit required a backward integration in time to 
obtain a solution. This, in turn, required a weight correction to establish 
the correct initial T/W q when the solution converged. These difficulties 

were avoided in the analytic solution by solving the retrograde problem in 
the forward sense, and only the previously mentioned difficulties were en- 
countered. An open-loop solution was initiated at a time-to-go of approxi- 
mately 10 seconds. 
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Altitude, h, ft/sec 
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( c) Altitude . 
Figure 5.- Continued. 
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(d) Pitch angle. 
Figure 5 • - Concluded . 


TABLE II.- LANDING FROM LUNAR ORBIT 
t/W q = 0.4; I = 315 sec; y Q = 50, 0°° v 0 = 5674.5; r Q = 0 


sec 

Q q , deg 

0 1# d.eg 

V , ft/sec . . . . 

y x > « 

v , ft/sec . . . . 

r- L> d-eg 

X t ) ^ 


Optimum, integrated 

290.02 

181.87 

163.23 

4,658.9 

15,000.0 

1,113.2 

-11.13 

1,019,285 


Analytic, guided 

290.92 

182.93 
164.28 

4,677.0 
14,998.8 
1,114.3 
-11.07 
1, 022, 944 



Figures 6(a) to 6(d) are the time-history comparisons for a guided in- 
sertion into a circular earth parking orbit. Table III presents the numerical 
results at the boundary. In this trajectory, velocity to be gained rather 
than time-to-go was the terminating criteria. This necessitated thrusting for 
an additional 0.13 second to meet the desired end conditions. This shift in 
the terminal boundary condition is readily apparent in figure 6(d) where it is 
observed that the pitch angle rapidly changes value. 
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26.0 x TO I * 3 


I ime, sec 


(a) Velocity. 


Figure 6.- Optimum insertion into a circular earth orbit 
T/W 0 = 0 . 751 ; I = 420 sec. 
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0 40 80 120 160 200 240 280 320 360 400 440 


Time, sec 

(d) Pitch angle. 

Figure 6.- Concluded. 

TABLE III.- INSERTION INTO A CIRCULAR EARTH ORBIT 


T/w 0 = 0.751; I = 420 sec; y 0 = 181,500 ft; v Q = 8,750 ft /sec; r Q 



Optimum, integrated 

Analytic, guided 

t^, sec 

417.15 

419.07 

9 0 , deg 

33.73 

32.21 

9 X , deg 

4.04 

1.38 

V , ft/sec .... 

18,537 

18, 733 

ft 

600, 932 

600, 930 

ft/sec .... 

25,570 

25,570 

r ± , deg 

0 . 0000 

-o.ooo4 

x, , ft 

5,929,025 

5,960, 834 




Trajectory Scanner 


Figure 7 indicates a parametric study in T/W q of landings from 

positions off pericynthion. The state variables and performance data for 

changes in the true anomaly angle are presented in table IV. A point of 

interest in this curve is the small penalty loss in performance for a wide 

variation in the true anomaly angle T| . For instance, in the case of 

T/w q = 0.4, a penalty loss of AV c = 10 ft/sec results for a range of T| 

of approximately 50° . For reasonable values of T/Wq, this small loss in 

AV indicates that lunar landings may be initiated from a wide range with 
c 

only a . small penalty in performance . 
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True anomaly, V , deg 

Figure 7 .- Landings from positions off pericynthion. I = 315*0 sec 



TABLE IV.- LANDINGS FROM POSITIONS OFF PERICYNTHION 



= 315 sec; 


y 1 = 1,000 ft; v^ = 100 ft/sec; 


= 0 


*] 



Quasi-optimum guidance equations in two dimensions for which constant 
thrust was used have been analytically derived by the calculus of variations . 
The basis of the analysis is the analytical solution to the flat-body prob- 
lem which is extended to simulate a round body by the addition of power series 
in time. Numerical solutions to the resulting transcendental equations have 
been generated by using an iterative convergence technique. 

Guided trajectories were examined for a landing on and launch from the 
lunar surface and for an insertion into earth orbit. The resulting curves for 
the control and state variables agree closely with the true optimum and pro- 
duce solutions which are more linear than the exact curves. The trajectories 
show losses in characteristic velocity of 0.564 percent for the lunar launch, 
O.389 percent for the lunar landing, and 1.057 percent for the insertion into 
earth orbit. This last trajectory had a terminating criteria on velocity 
rather than on time-to-go which had the effect of increasing the thrusting arc. 


Manned Spacecraft Center, 

National Aeronautics and Space Administration, 
Houston, Texas, February 25, 1964 



APPENDIX A 


OPTIMUM SOLUTION WITH FUNCTIONS OF THE STATE VARIABLES PRESCRIBED AT A 
VARIABLE TERMINAL TIME, INCLUDING MINIMUM TIME PROCESSES 


The material contained in this appendix was taken from a lecture series on 
optimization of dynamic systems at Harvard University in 1963. 

Given a system of first-order differential equations of the form 

q = f [q(t),ett),t] q(t 0 ) given t Q < t < t ± (Al) 


where <q(t) and 0 *(t) are the state and control vectors, respectively. It 
is desired to minimize some performance index of the Mayer form 


I 



(A 2 ) 


Certain boundary conditions are desired of the form 


i? 



= 0 


(A3) 


The performance index and the boundary conditions adjoined to the constraint 
relationships (eq. (Al)) with Lagrange multipliers produce the integral 


I = I G + u 


i? 


t=t. 




q(t),ett),t 


- q(t)\ dt 


( A4) 


where v* are constant Lagrange multipliers and the superscript T refers to 
the transpose matrix. 


The modified Hamiltonian is defined as 


H = A* T f [ q(t),ett),t] (A5) 

The differential of (A 4 ) taking into account differential changes in the ter- 
minal time t^ is 
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1 


t t /, sg -»t .. , 

dI = < l St + u 5t dt + 


Bg -»t Be? 

Sq(t) + V Bqjft) 


d?(t) 


t-fc. 


5§8qXt) +^8e > (t) - ?? T Bq^t) 


0 


dt 


(A6) 


Integrating by parts and using 

a ?(‘i) * s «(h) + ?( t i) dt i 

then. 


dl = 


BG , -*T Si? , -4/.x 

5t + u St + A q(t) 


dt ;L + 


BG , -» T Be? =r>T 

Sqjt) U Sift) 


dqtt) 


(AT) 


i. 


Sq(t) 


.*1 r ' 


t=0 


BH T \ R -» , BH erf 

Si + A ) 5<1 + ^- sr 


;t-t 1 


dt 


(AS) 


The Lagrange multipliers are chosen such that the coefficients of Sq., 
dq(t ^) , and dt^ vanish if t^ is not defined. 


^T 


T Bf 


* ss 


(A9) 


A 


TU \ /SG ->T BeT 

k 5J + u Si 


'(l) - 


(A10) 


t=t. 


Bg -»t BE? , _ 

St +v St +H 


= 0 


t=t_ 


(All) 


As a result of this choice of A(t) equation (A8) is reduced to 


dl 


■( 


=r*T r ->\ 

A 5 U=o + 


S’ & ldt 


(A12) 
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For a stationary value of I 


8h t 3f 

5§- = x 5e = 0 


(A13) 


For minimum fuel consumption or under the assumption of constant thrust for 
minimum time t i, = t^ such that 

'da' 


5tJ = 1 

't=t ± 




CL =0 

t-ti j 

For the specific case under consideration, the functions f q*(t) ,9*(t) , t 

Q ' — -J 


(AliO 


are 


x = u 


y = v 


u = ^ cos 0 


■ 'p • „ 

V = ~ sm 9 - g n 

(j. 0 


(A15) 

(A16) 

(A17) 

(Al8) 


and the boundary value constraints j? £q.(t^,tj 

■(h) - y ! - 0 


are 


y 


i\) -i 


u = 0 


V 


(h) - -i 


v, = 0 


(A19) 

(A20) 

(A21) 


where allowed to be free. 
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Applying equation (A9), the Euler- Lagrange equations are 


\ = ° 


a 2 = 0 


" “ \ 


\ " " A 2 


(A22) 

(A23) 

(A24) 

(A25) 


Applying equation (A13) to obtain the optimum control laws 


tan 0 = A^^y A ^ 


(A26) 


Integrating equations (A22) to (A25), and using the results in equa- 
tion (A 26 ) 


Therefore^ 


Applying equation (A10), 



Cit 


c 2 t 


C 2 t 


Cit 


Ai(ti) = 0 

A 2 ( t 1 ) = u 2 
^(h) = 
\{\) = \ 


(A27) 

(A28) 

(A29) 

(A30) 


(A31) 

(A32) 

(A33) 

(A34) 

(A35) 
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Ill' 


From equations (AJ2) and (A27) it is observed that the constant = 0. 


The optimum control law reduces to 


tan 8 = 


K - ^ 


One more boundary condition results from this analysis. Applying 
equation (All) along with the condition (Al4) 


1 + A 2 y + A^u + A^v 



APPENDIX B 


INTEGRATION OF A SET OF LINEARIZED DYNAMIC EQUATIONS IN WHICH A CONSTANT 
THRUST AND A LINEAR TANGENT LAW OF STEERING ARE ASSUMED 

Given a system of equations described by 


x = u 


y = v 


6 

u = ^ cos 0 - / h _ t 

M* Li n +l 

n=0 


n 


(Bl) 

(B2) 

(B3) 


v = sin 9 + / a t n 

M- Z__, n+1 

n=0 


(B4) 


tan 0 = 


% - c 2 t 


= Z 


(B5) 


along with the subsidiary equations 


= 1 + )j.t 


qt n 


<x„ = 


2 tan 0 - tan 9^ 


= 1 _ a 2 tan 0 Q 


k = cx^ + a. 


2 2 
'2 


(b6) 

(B7) 

(b8) 

(B9) 



sin 0 = Kz/( 1 + z 


cos 9 


(BIO) 
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K = 1 
K = - 1 


(-at/2 < 0 < Jt/ 2 )] 
(at/2 < 0 < 3at/2) 


(Bll) 


it is required to integrate equations (Bl) to (b 4) to obtain three simultaneous 
algebraic equations. Rewriting equation (B4), 


u 

ydv - |o> Ivx t n dt = cp sin 0 dt 


n=0 


(B12) 


integrating 


Z ,n+l _ 

v -1 a .nt 


v> <m- v - p-y 


n+l 


n + l 


n=0 




n + 2 


= a, (cp sec 0 - A) (B13) 


where 


A = - v o) + 'P sec e o 


(B14) 


Equation (B13) is of the form 


y + P(t)y = Q(t) 
-2 


(B15) 


and has an integrating factor of y. 

Dividing equation (B13) by y, 2 and noting the exact differential yields 


• d 

^ dT 


y/d - 


1 . n+2 

a t.. 
n+l 1 

y(n+l) (n+2) 


n=0 


CL^ 


= -+■ (cp sec 0 - A) 


(Bl6) 


The left side of the equation integrated between the bounds at t = 0 and 


t = 


t^ is 


t , n 

< V Vah 

y iri ' y 0 ' L Hitn+l) 
n=0 


4-2 


p,- L (n+l) (n+2) 
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The first term on the right is 


a^A 




- a 2 A 


■o 


at _ 

2 • J '2 

l-i M 1 


omAt_ 


= _ 2 1 


The second term is 


a 2 T 


see e at r tan 0 1 ( 1 + it!l 

^ " i in,„ ^ + “ 2Z )‘ 


dz 


_ Kc^ 

M- 


log^ 


z + 1 + z 


2 V 


1/2 


- a 


(les 2 ) 1/2 

2 a 3 _ + a 2 2 


i -l/2 , 

- ^ k log e . 


k' 


“1 + a 2 


r- ■* 

.ri 

\ 

) - 2a 2 ^l + z 2 ] 

2 


z 

1 

/ 


J 

/ 


, tan 0. 


tan 0 


sec 0 n 


„ , /tan 0 + K sec 0 \ / 

= ^<log (- 1 • ) + Ka sec 0 - 

^ | e \tan 0 Q + K sec Q Q J 2\ 0 p, 


- c^k -1 / 2 log 0 


1/2 


- cc^ - Kk 7 sec 9^ 

7 i / 2 


^k - - a g Kk 1 / 2 sec 9 q) 

where the following identities are true: 

p,(t = 0) = 1 

a i * a 2 tan 6 1 = ^1 


0 


a i + a 2 tan 6 0 = 1 





Ti. • icegral of equation (Bl6) is 


l . n+2 

v -1 a .,t. 


V n+l“l Kp J 
A -(n+lHn^) f < los . 


tan 0 1 + K sec 0^ 


el tan 0 Q + K sec 0 O 


+ Ka o( sec - sec 


- a^" 1 / 2 log e 


k - sec 9 1 

( k C, a i " a- 2 Kk 1/2 sec 0 O )^ 


A a_t 


- 2 t i/M 


(B17) 


Integrating equation (B3) 


u l - u 0 = 9 


cos 0 dt 




h _ t dt 
n+1 


0 n=0 


•Kk“ 1//2 cpa 2 


Using equation (Bl8) in (B17) 


f 1 „ + n+2 

•2 V a n+l t l 

q qqy 0 - y ± + ( n +l)'(~2T 

L n=0 

/ sec e 1 

+ KoJsec 0 Q - — — 


k - ot-qM-q - OgKk / sec ®q 


k - - agKk sec 0 o Jqq 


(B18) 


/ tan 0^ + K sec 0^ 
[j,A a 2 t ± = - Kipqq lo S e ^ tan 0 Q + K sec 0 Q 


cp a 2 \ 1 


u 0 + *1 


(B19) 


where 


♦1 = 


l , , n+l 

-1 h , n t n 

S n+l 1 

\ n + l 


(B20) 


Equations (Bl8), (B19), and (B13) evaluated at the Boundary t = constitute 
the necessary three algebraic equations in the three unknowns t^, 0 q, and 9q* 
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APPENDIX C 


FIRST-ORDER PERTURBATION TECHNIQUE FOR THE SOLUTION 
OF TRANSCENDENTAL EQUATIONS 


Given a system of algebraic equations defined by 

f i( 6 j)"° 


(Cl) 


where the 9 . are a set of unknown variables, it is required to determine the 
unique set of 0^ such that equation (Cl) is true. The first-order perturba- 
tion equation is: 

c)F. 

"i ~ ( F i)m + 1 - ( F i)m - 557 < C2 > 

d 

where m defines the order of correction • 

Define the matrix 


It is desired that 


dF. 

^ = a ±j = SeT 


( F i)m + 1 = 0 


(C3) 


(04) 


Using equations (C3) and (c4) in (C2) and solving for the correction term 


ae , -t - 1 • (F ± ) m 


(05) 


The functions F^0 ^ are defined as 

1 / 2 . 


-Kk 


F = 

1 cpa 


L ( U 1 " "o + *l) + l0g e 


l/2 ~ 

k - - ctgKk ' sec 0^ 


(k - - a 2 Kk 1 ^ 2 sec 


(06) 
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II 


Z i n+1 

a .,t. 


V Wl 

= - m. -wr 1 - ) — + i 


n=0 


1 + |J,t 


(n+1) 


1 n + 2 


+ CL, 


,(cp sec 0 1 - 


„ .2 
f 5 = H 


Z , n+2 

a , t. 


, V 1 a n+l 1 

x i y 0 " y l + (n+1) (n+2) 


n=0 


- |iA a 2 t 1 


+ Kcpp^ 


tan 0 + K sec 0 \ / sec 0. 

log It t — TT — ]+ Ka (sec 0. - 

e\tan Q n + K sec 0 n / 2\ 0 M., 


Kpa^ 

cpa. 


(”l - U 0 + *l) 


where 


A = ^ y o - v o) + * sec e o 


Z , n+1 

V Wi 

*1 * 1 ~+-l- 
n=0 


|at n 


= 


2 tan 0^ - tan 0^ 


= 1 - tan 0 q 


, 2 2 
k = a x + a g 


K = 1 


K = - 1 


(- it / 2 < Q fl / 2 ) 

( rt /2 < 0 < 3 rt/ 2 ) 




and the 0 described previously are t^, 0 Q , and 0^. 


A) (C7) 


(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 

(C14) 
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ii ii i 


ii i ii 



r 


a 


-s, 5a ,/p f 1/2 Sa 2 5k/B0 1 \ 

557 - “l 59^ ' V k S6C ®1 *“ 6 1 - ' 59 7 + a 2^~p72j SeC 


13 


1/2 

k - a^ - a^K k ' sec 6 


Sk 3a l L 1/2 9a 2 - Sk/Se i 




- Kk 


v + CL_ 77 0 1 s e G 0^ 

^®1 2 2K k 1 / 2 / 0 


172 

k - - a^K k ' sec 0 Q 


Kuk -1 ^ 2 / . \^ a 2 5k , ^ a 2 


a ) 


5a^. 


5a •/• 


a 21 = ^ sec 9 1 " A ]5i~ “ a 2^ - ^ V 1 “ M 'l a j) 


da, 


22 


/ A '2 5A 

= (cp Bee 0 1 - A)^- - a 2 


da 


dA 


23 


/ \ — 2 / a A 

= ^cp sec 0 1 - AJ^~ + ajcp sec e ± tan 0 1 - ^ 


a 31 = ^ 


5a^ 


•/. \ A ^ /. 2 aA 

n(w 0 + 0 l) - a 2 A ' *1^* 5t^ + “2St7 

/ tan 0 n + K sec 0 


+ cpn 


|_ia. 


1 /.. 


K loe 1 — " — ■ — r— • ] + a„ sec 0^. - — ~fu n - u_ + i|r, 

S e V tan 0 q + K sec ®q / 2 ® ^ a 2 \ 10 1 


sec 0A5a 2 ^ 


+ ^iS\ sec e o - ~vr)^\ " 


°i +3 - 


( U 1 - u 0 + *1 


(C17) 

(C18) 

(C19) 

(C20) 

(C21) 
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J>A . . aa 2 


a j2 - - ut il a 2 5e7 " A 557 ) + w i 


0 / 


sec 0^ - 


sec 0-jN da^ 


0" n x / ^9 


0 


da. da. 


- a, sec 0^ - 


- 0 2 
9 a 2 


fu x - u Q + f^k - a. 




ZA . H 
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= - a 2 ^07 + A 5©7 + 


sec 0 \ da^ 

sec 9~ - \ S T ," 

o n x pi 


CL \ 

+ [ 1 - — tan 0., i sec 0 n + 
1 1 


da^ 


cp a 2 


2 ( U 1 ' U 0 + +l) ^ 


where 


a i = 


l , n+l 

a t 
n+l 1 

n + l 


n=0 


l . n+2 

a t, 
n+l 1 


CT 2 /_, (n+l) (n+2) 

n=0 


<T = / O' t_. 

3 n+l 1 

n=0 


n 


% Z_. h n+l t l 

n=0 


n 


So- 


5tT = V 1 


* 1-0 


(C22) 


(C23) 


(C24) 


(C25) 


(C26) 


(C27) 


(C28) 



(C29) 


Sa- 2 


a 


'2 

M-t-. 


2 

sec 9. 



2 

a 2 2 

sec 0., 

• _L 

[it ± 


(C50) 


SA kx 

5tT = “ 


a*. 




(C51) 


Sa 

tL_ ( v _ ny '\-sr— — • + cp sec 0 „ tan 0 (C 32 ) 

^ - 2 V V 0 Wo)XT 0 ^00 


Sa 

30^ 




(C33) 


Sen Sa 2 

= “ tan 9 o 


Sot, Sen 2 

^ = “ tan 9 0 ^ “ a 2 S6C 9 0 


So, Sa 2 

S5^ ‘ * tan 9 o 557 


S°°2 

zq - 2 ( a 2 ' “i tan e o) s*I 


(C3^) 

(C35) 

(036) 

(C37) 



Sk 







2a l “2 


2 

sec 


(C58) 


(C39) 


Of the three variables concerned, t n , 0^, and 0^, only the latter two re- 
quire initial guesses to initiate the solution. A minimum value of time-to-go, 
t , may be computed from the following relationship. The minimum characteris- 
tic velocity is: 

M . -l AV l = cloe efit) 

' 'mxn 1J 


where AV is the change in the total velocity and c is the effective exhaust 
velocity. Now, 


= 1 + nq 


(o4i) 


Eliminating |j.rb from equations (C^O) and (C4l) and solving for t^. 



min 



(042) 


One more point requires definition. It is desired that the functions 
^ be driven to 0 by means of the corrections dq^ as stated in equa- 
tion (Cb) . Since for computers this method is impractical, tolerances (or 
errors tolerated) should be defined. Only one tolerance requires definition 
since the variables are related by equations (Oil). When the tolerance on 
time-to-go t is defined as e . , the following tolerances are derived: 


e 



(043) 
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( 044 ) 
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